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it 0 - DECAY Y-RAYS FROM THE GALAXY 
AND THE INTERSTELLAR GAS CONTENT 

F. W. Stecker 

Theoretical Studies Branch 
Laboratory for Space Physics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

Abstract: Knowledge of the total y-ray production rate 

per H atom from the decay of tt^ mesons produced 

in interstellar cosmic-ray interactions is 

essential for determining the possible amount of 

interstellar H 2 . This production rate is 

recalculated here using the latest accelerator 

data on production in pp-interactions up to 

^300 GeV. A simple but accurate approximation 

used here resolves the past disagreement over the 

magnitude of this rate. An upper limit is 

- 2 5 - 1 

obtained of (1 . 51+0 . 23) xlO s consistent with 

. . - 2 5 - 1 

the observed upper limit of 1.6x10 s 

Knowledge of the total Y-ray production rate per hydrogen 
atom from interstellar ir G -decay is essential for determining 
whether molecular hydrogen may exist in significant quantities 
in interstellar space. Recently, Kraushaar, et al. (1972) 
have obtained an observed galactic Y-ray flux measurement 

_ 2 5 

corresponding to a production rate of 1.6x10 photons/s 
above lOOMeV for each interstellar hydrogen atom observed 
in 21cm radio emission. At the same time, recent theoretical 
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calculations have given diverse results on the Y-ray 
production rate per H atom. These values of 1.3xl0 _25 s _1 
(Stecker 1970), 1.8x10 2 5 s ^Cavallo and Gould 1971) and 
3.2x10 2 5 s 1 (Levy and Goldsmith 1972) are enough in disagree- 
ment with each other to influence qualitative conclusions on 
the interstellar medium. If we compare the value given by 
Stecker (1970) with the observations we conclude that there 
may be some contribution to the observed y-ray intensity 
from "cool" H and H 2 not observed in 21cm radio emission. 
(Stecker 1969, Stecher and Stecker 1970). The rate given 
by Cavallo and Gould (1971) would appear to allow for no 
significant contribution of cool atomic and molecular hydrogen 
to the galactic y-ray production rate and, even so, would 
appear to be a bit high. The rate given by Levy and Goldsmith 
(1972) is so high as to throw the observation of Kraushaar, 
et al. (1972) into question. 

The three calculations referred to above involved elaborate 
numerical integrations to determine first the differential 
Y-ray spectrum and cannot be easily checked by the reader. 
Fortunately however, if the calculation is limited to the 
total y-ray intensity, a simple treatment can be given which 
the reader may check for himself. In this manner, a reliable 
figure for the y-ray production rate may be arrived at. 

Figure 1 shows an up-to-date summary of the accelerator 
data on the total cross section (a) times multiplicity ( C) 
for neutral pion production in p-p interactions for energies 
up to % 300GeV shown as a function of kinetic energy (T) as 
given in the references. These data are well approximated 
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by the broken power-law 


% o (T) (T) 




r 


10 - 2 5 T 7 • 64 cm 2 

8.4 x 10“ 27 T° * 5 3 cm 2 


0 . 4^T^0 . 7GeV 


T^-0 . 7GeV 


as the reader can verify from the figure. Taking the cosmic- 
ray spectrum I(T) = 0.15 T~ 2 * 2 cm" 2 s _1 Sr“ 1 Gev“ 1 used by 
Levy and Goldsmith (1972) , the total y-ray production 
rate from p-p interactions is given by 


q yH = 8 tt J dT I(T)a i0 (T)^o(T) 

0.7 

= 3.77xl0~ 25 / T 5 ’ 4 4 dT + 3.17 x 10 2 6 / T~ 1,67 dT 

0.4 0.7 

= 0.66 x 10~ 2 5 s~ 1 


Adding in the contribution from p-a, a-p and a- a 
interactions in the galaxy brings the total production rate 
per hydrogen atom up to '\JL0 -25 s~ 1 , a value in good agreement 
with Stecker (1970) (which used a somewhat higher cosmic-ray 
intensity) and with previous work summarized by Stecker (1971) . 
There is, of course, some uncertainty in the assumption of 
the true "demodulated" galactic cosmic-ray spectrum as 
distinguised from that observed at the earth. However, 
using the upper-limit to the demodulated cosmic-ray spectrum 
given by Comstock, et al.(1972), an upper limit on the y-ray 
production rate is obtained of (1 . 51+0 . 23) xlO 2 5 s 1 with 
the error bracket reflecting the experimental error in the? 
accelerator data on a£. The above value is consistent with 
the value of 1.6xl0 -25 s -1 given by Kraushaar, et al.(1972), 
which also represents an upper limit since it does not take 
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account of the additional contribution from cool H and H 2 
which may be adding to the observed flux. 

It is hoped that the straightforward calculation pre- 
sented here will at last resolve the discrepancy existing 
between various theoretical calculations and finally result 
in a consistent picture of y-ray production in interstellar 
cosmic-ray interactions. It should, of course, be noted 
that whatever the shape of the differential y-ray production 
spectrum (and here there is more uncertainty because it 
depends on the dynamics of pion production) , the normaliza- 
tion to the total production rate has to be consistent with 
data on the total cross section and multiplicity. 
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